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RELATIVE TRUE AMPLITUDE MIGRATION 

BACKGROUND OF THE INVENTION: 

Conventional "true-amplitude migration" is a form of migration of the Kirchhoff 
type using a weighted diffraction stack. This technique has become increasingly popular 
as a seismic tool for quantitative discrimination of lithology and fluids in hydrocarbon 
reservoirs. However, the term "true amplitude migration" is to some extent 
misdescriptive. From a physical viewpoint, valid amplitude information may only be 
preserved if the entire processing line is designed for that purpose. In most cases, it is 
not. Thus, it is virtually impossible to have a "real" true amplitude migration in practice, 
because many processing procedures affect the amplitude preservation. 

Furthermore, conventional true amplitude weight terms are complicated and to do 
not consider many necessary factors. In a smoothly varying velocity field, for example, 
the amplitude weight depends on the travel-time, amplitude, ray-tube spreading factor, 
and takeoff angles in the shot and receiver points. 

Even further, conventional approximations of the reflection travel-times often 
assume a small offset-to-depth ratio, and their accuracy decreases with increasing offset- 
to-depth ratio. Hence, they are not suitable for migration of long-offset reflection seismic 
data. 
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Finally, conventional Kirchhoff time migration for transversely isotropic media 
with a vertical symmetry axis (VTI media) is implemented using an offset-midpoint 
travel-time equation. However, travel-times as functions of offset and common midpoint 
(CMP) are given by the simple analytical double-square-root equation (DSR) designed 
for homogeneous media. This travel time is not well suited for inhomogeneous media. 

Thus, there is a long felt need for a relative true amplitude migration process 
which preserves relative amplitude, considers the appropriate factors, is accurate for long- 
offset data, and is suitable for vertically transversely isotropic (VTI) media. The present 
invention addresses the above mentioned concerns. 
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SUMMARY OF THE INVENTION: 

In one example embodiment of the present invention, a method is provided for 
migrating seismic data. The method comprises determining a travel time; determining a 
ray-bending correction; determining an amplitude correction; and outputting data based 
on the travel time, the ray-bending correction, and the amplitude correction. 

In a further embodiment, a system is provided for migrating seismic data. The 
system comprises means for determining a travel time, means for determining a ray- 
bending correction, means for determining an amplitude correction, and means for 
outputting data based on the travel time, the ray-bending correction, and the amplitude 
correction. 

In an even further embodiment of the present invention, seismic data is provided. 
The seismic data is produced by the method comprising determining a travel time; 
determining a ray-bending correction; determining an amplitude correction; and 
outputting data based on the travel time, the ray-bending correction, and the amplitude 
correction. 
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BRIEF DESCRIPTION OF THE FIGURES: 

FIGURE 1 is a diagram of a horizontally layered medium of an example 
embodiment of the present invention. 

FIGURE 2 is a diagram of the geometry of an example embodiment of the 
present invention. 

FIGURE 3 is a diagram of a seismic gather using a conventional migration 
(straight ray) method. 

FIGURE 4 is a diagram of a seismic gather using an example embodiment of the 
present invention. 

FIGURE 5 is a diagram of a seismic gather using a conventional migration 
(straight ray) method. 

FIGURE 6 is a diagram of the outputs of an example embodiment of the present 
invention which accounts for transversely isotropic media (VTI). 

FIGURE 7 is a diagram comparing amplitudes of an example embodiment of the 
present invention and a conventional migration technique. 

FIGURE 8 is a diagram of a model. Densities are displayed although velocities 
and densities were used for the seismic modeling. 

FIGURE 9 is a diagram of the outputs of a conventional migration technique 
using the straight ray approximation of model data corresponding to the model from 
Figure 8. 

FIGURE 10 is a diagram of the outputs of an example embodiment of the present 
invention. 
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FIGURE 11 is a flow chart illustrating a job flow of one example embodiment of 
the present invention. 
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DETAILED DESCRIPTION OF EXAMPLE EMBODIMENTS OF THE 
PRESENT INVENTION: 

In one example embodiment of the present invention, a method of Kirchhoff 
migration is presented that includes ray-bending and relative true amplitude in 
transversely isotropic with a vertical symmetry axis (VTI media). The new method is 
called "relative true amplitude migration" (TAPSTM). In one embodiment, the method 
uses a Kirchhoff operator and determines the kinematics or the shape of the operator from 
travel-times computed assuming V(z) and VTI media. In a further embodiment, to 
include ray bending, the earth is assumed to consist of horizontal layers so that the 
velocity profile used to calculate the travel times is simply a function of depth or time. 
Because ray-bending and VTI are considered, the takeoff angles for shots and receivers 
are calculated accurately compared to those obtained assuming straight rays. Thus, in 
further embodiments, the amplitude correction terms of the operator are determined using 
the takeoff angles for preserving the relative true amplitude. The method is robust in the 
presence of an irregular geometry and therefore lends itself well to application on most 
modern seismic marine, OBS, and land geometries. 

In a further embodiment, the method has the following main steps: 1. A time 
table is constructed to calculate the output imaging time. 2. The input is duplicated in the 
time domain for anti-alias filtering and resample. 3. The velocity, VTI parameter and 
mute tables are set up for each CDP and distributed uniformly among the processors. 4. 
The migration starts. The first processor reads a trace (record) and broadcasts the trace 
(record) to all processors. Then, all processors begin imaging the data. 5. The migration 
operators are built by calculating the high order (ray-bending effect) and VTI travel-times 
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and amplitude correction terms. 6. TAPSTM completes the migration and outputs data in 
gather or non-gather modes. 

Ray-bending for a large offset-to-depth ratio 

In a further embodiment, the method further comprises determining a ray bending 
effect for large offset to depth ratios. Those of ordinary skill will recognize that three- 
dimensional seismic reflection data is acquired on an increasingly routine basis. These 
large volumes of data offer the potential for very high-resolution 3D images of subsurface 
geology and subsequent estimation of the earth's physical properties (Lumley and 
Beydoun, 1992). In seismic acquisition, energy from all sources is spread from one 
scatter-point to all receivers. Consequently, all input traces can contain energy from the 
scattering point. When the input traces have a finite recording time, the scattered energy 
is restricted to traces within the prestack migration aperture of the scattering point. The 
objective of prestack migration is to gather this energy and sum it back to the scattpr- 
point location. The key points are how to obtain a accurate travel-time for having better 
imaging and preserving the relative true amplitude. 

The double square root equation to compute the travel-times is fundamental in 
conventional Kirchhoff migration. In these conventional migrations the velocity is 
assumed to be constant, then straight ray is used for travel-time computations at an 
imaging point. Thus, the velocity change in vertical and lateral directions should not be 
too large. The conventional double square root equation is exact in the sense that there 
are no error terms dependent on dip angle and offset. However, this conventional double 
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square root equation assumes straight ray paths through the earth and does not faithfully 
correct for ray path bending associated with vertical changes in velocity. This is 
particularly erroneous for wide-angle ray paths such as those at steep dips or long offsets. 
Ray tracing travel-times include these effects and are more accurate. In a further 
embodiment of the present invention, the double square root formulation is improved by 
including higher order move-out terms. In principle, the move-out equation for V{z) 
consists of an infinite series of higher order terms. The number of terms needed in order 
to be accurate depends upon the vertical velocity profile and the propagation angle. The 
ray tracing approach removes these issues by accounting for the bending accurately. 

Figure 1 illustrates a layered media. dk 9 vk is the thickness and velocity at Mi 
layer respectively. Thus, the down or up travel-time is expressed by: 



Tx = V— ~ minimum (1) 

Pl=d 2 k+ x 2 k (2) 

n 

* = E** ( 3 ) 



Jt=l 



Where x = x r or x = x* 9 which is the distance between the receiver and imaging point, 
and the shot and imaging point, respectively. The minimum requirement in equation (1) 
under the constraint of equation (3) leads to: 

dTx x. 



A (4) 



dxk p k v k 

where lamda is a Lagrange multipler. From equation (2) and equation (4) we get 
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2 ^y>\d\ 

x *=^kr ( 5 ) 



where 

2 2 
X k _ X k ^2 

P>1 (d 2 k +x 2 k) v 2 k ~ 

Thus, we have 

v.d 



x ■ 



*Z , ' 2 2 (6) 



T-t-jLIZ- (7) 

Equations (6) and (7) are the parametric forms of the time-distance relationship 
with Lamda as parameter as given in Slotnick. 1 As Lamda varies from 0 to 7/v, where 
v=max(vl } v2, vw), x will vary from 0 to infinite. We have 
sin 0. 

JL = k - (8) 

Let 

7/ =Cj +c 2 x 2 +c 3 x 4 +c 4 jc 6 +c 5 x 8 +c 6 x 10 + •• (10) 
Where we have 

_ 2 _ fl i _a2- a i a 3 _ 2aj«3 -aja 2 a 4 -«2 a 3 
Cj — a 1 c 2 — c 3 — - - c 4 — ■ 



a, 



a x 4a 2 
= 2 Z^ «2=2JM a 3 =2±vld k a 4 =2±vld k 

k=\ V k k=\ k=l k=l 



1 Incorporated herein by reference. 
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In still a further embodiment, to obtain a travel-time approximation with good 
accuracy at large offsets, a new series which is based on a large offset-to-depth ratio 
rather than a small one is derived: 



+ c 2 x 2 + c 3 x 4 + c 4 x 6 + c 5 x s + c 6 x 10 + • • 



- \{c x + c 2 x 2 + c 3 x 4 )(l + , ) 



= -yj(c i +C 2 X 2 +C 3 X 4 ) 



2 4 

^ 1 c,x 6 +c,x s +c fi x 10 +•• ^ 

1 + — - 5 — — - — - — +•• 

2 C-y ~\~ C 2 X ~h C^X J 



(11) 



c 4 x 



1 + -CC- 2 4 

y 2 ~h C 2 X H~ C^X j 



CC is a constant which depends on the new series. Thus, the total down-up 
travel time is determined; 



T = T g +T S 



go 



( 1 ^ 
2 T 



sO 



( 1 

2 T 



(12) 



In equation (12), the fourth and sixth order series correspond to the large offset 
asymptotic travel-time. This travel-time has a higher accuracy than conventional travel 
times using the double-square-root equation. In addition, equation (12) is derived from 
the general sixth order series and it should be more accurate. Thus, in a further 
embodiment, equation (12) is used to calculate the travel-time for the large offset-to- 
depth ratio. Later, we will demonstrate how the travel-time calculated with equation (12) 
is close to the true travel-time. 



Travel-time in transversely isotropic media 
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In still a further embodiment of the present invention, a Kirchhoff time migration 
for transversely isotropic media with a vertical symmetry axis (VTI media) is provided 
using an offset-midpoint travel-time equation. The derivation of such an equation for 
VTI media requires approximations that pertain to high frequency and weak anisotropy. 
However, the resultant offset-midpoint travel-time equation for VTI media is highly 
accurate for even strong anisotropy. In one embodiment, the travel-time equation is 
dependent on two parameters: migration velocity and the anisotropy parameter (eta). 

In homogeneous VTI media, it has been demonstrated that P-wave velocity V P o 
and travel-time are practically independent of S-wave velocity, even for strong 
anisotropy. This implies that, in a further embodiment, for the migration purposes, the 
travel-time is expressed, 



where 8 and 8 are anisotropic parameters. 

The travel-time calculated with equation (13) is accurate. It is also stable for a 
large offset-to-depth ratio. 




(13) 



where, 



7 = 



s-6 
1 + 2S' 
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Thus, in a further embodiment, the total down-up travel-time for a VTI media is, 



T = X + - c(vu)x: + r 0 2 + i - c(vii)x: (14) 



where 

2 7 



C(H7) = 



The offset-midpoint travel-time formula, equation (14), for VTI media is similar 
to the double square root equation for isotropic media but as a function of the source and 
receiver ray parameters. In addition, the efficiency in travel-time calculation using 
equation (14) allows its practical use in parameter estimation. 



Travel-time with ray-bending and VTI media for a large offset-to-depth ratio 

Equation (14) is a formula to calculate the travel-time for a VTI media at small 
offset-to-depth ratio. In conventional methods, we use ray-tracing to calculate the travel 
time for a large offset-to-depth ratio for VTI media. This is a depth migration procedure 
and it is very time consuming. Here, in a further embodiment, we develop a new method 
to take VTI into account for time migration without ray-tracing. The travel-time 
difference between straight ray (with the double square root) and VTI is calculated, 

V rms \ v rms V rms \ rms 

Because the traveltime caused by VTI media is much shorter than that of the 
straight ray general case, in still a further embodiment, we make an approximation for 
simplifying code implementations, 
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AT ajrn)x A t ^ c{VTi) x A r 



(15) 



2jr 0 2 +4- 



iat: + 



vl 



Equation (15) is the travel-time difference between the straight ray assumption 
without VTI and with VTL In a further embodiment, combining equations (15) and (12), 
we calculate the travel-time to include both ray-bending and VTL 



T=T +T + AT 

5 -» 



= T 



( 1 

i + ±cc C4 * 6 



v 



T 1 



sO 



J 



i+-cc^ 

2 Tl j 



+ AT(VTI) 



(16) 



The difference AT is due to VTI, therefore, we combine equations (12) and (15) 
for VTI and large offset-to-depth ratio (equation 16), Equation (16) takes both ray- 
bending and VTI effects into account for large offset-to-depth ratios. 



Weighting function for the relative true amplitude migration 

Knowledge of elastic parameter (compressional and shear velocities and density) 
contrasts within the earth can be correlated to lithology and fluid changes. Elastic 
parameter contrasts manifest themselves on seismic records as a function of angle- 
dependent reflectivity. Interpretation of angle-dependent reflectivity, or amplitude 
variation with offset (AVO), on immigrated records is often hindered by the geometrical 
spreading loss, as well as other factors. Various conventional migration algorithms 
involving weighted diffraction stacks proposed by geophysicists are based on Born 2 or 



2 Incorporated herein by reference. 
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Kirchhoff 3 approximations. In a further embodiment of the present invention, the 
principal issue in the attempt to recover angle-dependent reflection coefficients becomes 
the removal of the geometrical spreading factor of the primary reflections. In one 
embodiment, the weight function that achieves this aim is independent of the unknown 
reflector and correctly accounts for the recovery of the source pulse in the migrated image 
irrespective of the source-receiver configurations employed and the caustics occurring in 
the wave-field. Thus, one principal aim of the relative true-amplitude migration consists 
of removal of the geometrical spreading factor from seismic primary reflections without 
any knowledge of the search for subsurface reflector. In a further embodiment, the 
weight formula derived by Schleicher et al. 4 is expressed as, 



Vcos^cos a g | det(rj Nsr + T^Nsr) 



(17) 



v, V| det(A^)lVl det(AW)| 



Now let us discuss the assumptions in equation (17). 



1 



The source pulse is reproducible. 



2 



The source and receivers fall on a planar surface, and reflection points lie on 



smoothly curved interfaces. 



3. 



The wave propagation is described by zero-order ray theory. 



4. 



The ray amplitude varies slowly with offset. 



5. 



The weight function does not vanish at the stationary point. 



Now let us explain all quantities that appear in equation (17), 



1. 



Angles as and a% are the takeoff and emergence angles of the central ray. 



3 Incorporated herein by reference. 
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2. Parameter v s is the P-wave velocity at source point. 

3. Ys and To denote certain constant 2x2 matrices describing the source and 
geophone locations with respect to the position vector ^. They depend on the 
measurement configuration. 

4. Nsr and Ngr are second-order mixed-derivative Hessian matrices of travel- 
times. 

It is important to note that the quantities appearing in equation (17) only depend 
on the two individual ray branches SR and RG and not on any reflector properties. In 
other words, they are computed for a complete ray SRG irrespective of whether or not 
there is a reflector at R. This has an important consequence. It means that equation (17) 
allows for a generalization to all points S and G specified by any coordinate pair (£1,^2) 
and any point M in the macro-velocity model. 

In a further embodiment, equation (17) is modified for coding. First let us 
suppose that the ray branches SR and RG are straight, thus, 

Vcosc&cos a z 



Vs 



T T 

8 ! *s 

s s 



4L 2 H 2 sm l (a) 



2rr4 



TT T V 

s g r • 



(18) 



Figure 2 illustrates the parameters that appear in equation (17). Tr is the imaging 
time, Ts is the travel-time from source point S to the reflector R, and T g is the travel-time 
from R to the receiver point G. Vr is the RMS velocity at the R. In still a further 
embodiment, equation (18) is derived based on the straight ray assumption. It might 
generate a large error in the weighting function calculation for a large offset H. As we 

4 Incorporated herein by reference. 
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know, the takeoff and emergence angles, as and a g , depend not only on the geometry, 
but also on the velocity model. Thus, in one embodiment, angles as and a% should not 
be computed as following, 

cos(#0 = — and coste) = — 
Ts T g 

Instead, using an alternate embodiment, we calculate the angles as and a g 

differently. We have shown that equation (16) gives an accurate travel-time, therefore, 

all quantities that appear in equation (18) are calculated using equation (16). Hence, we 

have 



r dTs^ 2 



cos(^) = 1- V s , and cos(^) = 

dxs) \ 



1- 



dTg I V] (19) 



Equation (17) is derived using the straight ray assumption, hence in some 
embodiments, the weighting function is in error. However, in alternate embodiments, if 
we compute the travel-times using equation (16) which accounts for ray-bending and 
optionally VTI, the calculation of the takeoff and emergence angles, as and Og using 
equation (19) should yield a small error compared to the error obtained using equation 
(18). 



Pre-stack Kirchhoff time migration with relative true amplitude and weak VTI 

media 

Kirchhoff prestack migration is based on the Kirchhoff boundary integral solution 
to the scalar wave equation. It is a mathematical statement of Huygens' principle in that 
the wavefield at any interior point in a volume can be reconstructed by a surface 
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integration of the boundary values of the wavefield, weighted by the Green's function 
impulse response to the wave operator. In some embodiments, Kirchhoff time migration 
is expressed by using shot (xs 9 y s ) 9 receiver (x g9 y g ), and imaging point (xi 9 yt) 
coordinates as 



Where I(T r9 x l9 y f ) is the image. T r is the image time. W(g 9 R) is a weighting function, 

D(T ^T s +T g ) is the time derivative of the input data, T is the total travel-time from 

shot point S to the receiver G going through the image point R. In a further embodiment, 
by substituting equations (16), (17) and (19) into (20), we have a formula for the 
Kirchhoff time migration with relative true amplitude and weak VTI media. 

Implementation 

Substituting OS z&dOR into equation (15), we have, 



Equation (20) is the kernel part of some embodiments of the time migration for 
VTI media with relative true amplitude including the ray-bending effect due to V(Z). In 
one example embodiment, FORTRAN is used to do the intensive computation, sample 
mapping from input to all of output targets, because FORTRAN is more efficient than C. 




(20) 



T= T +T 

■s g 




(21) 
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The kernel part, FORTRAN codes, takes over 85% of the total time consuming. C codes 
handle the I/O, parameter controls. 

In some embodiments, FORTRAN is used to do the intensive computation, 
sample mapping from input to all of output targets. In some embodiments, the kernel 
part, FORTRAN codes, takes over 85% of the total time consuming. C codes handle the 
I/O, parameter controls. 

Data Flow: 

In most embodiments of the present invention, the migration must be the last 
procedure in a job flow. It outputs stack or gather volume. Figure 11 illustrates a job 
flow of one example embodiment of the present invention. Unlike conventional 
techniques, example embodiments of the present invention have a logic to compensate for 
fold variations, the method has the capability of using multiple nodes, and will divide 
equally the user specified output target bins among the nodes so that each node will have 
every nth CDP as a target; where n is the node number. In most embodiments, the nodes 
must have sufficient memory to hold all of the target bins (including offsets when doing 
gathers). Insufficient memory will result in the job aborting immediately. The job will 
let the user know about how many processors would be required to have enough memory 
available. The compute time of the nodes does not scale linearly with the number of 
nodes. This means that increasing the number of nodes will decrease the elapsed time of 
the job, but the number of input traces processed per node hour will also decrease. 
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A one dimensional model was used for testing and comparisons between the 
conventional PSTM and this TAPSTM. Figure 3 shows the migration result by using 
conventional migration techniques. The gather output at the same CDP location is not 
flat. Figure 4 is the result of one example embodiment of the present invention. We see 
that the event at a CDP location is flat, which is what we expect. 

Meanwhile, we have tested a simple ID VTI model. Figures 5 and 6, 
respectively, show the results of example embodiments produced without and with VTI 
corrections. The travel-times are greatly effected by VTI. 

In order to test the amplitude preserving, we have generated two datasets by using 
the same model, the first one is without the spreading factor, the second one is with the 
spreading factor. We use an example embodiment, to migrate the second model data and 
output its imaging gather at the different CDP locations. The amplitudes of the example 
embodiment of the method for different offsets should be equivalent to that of the CDP 
gather after NMO to the first dataset. Figure 7 shows the amplitudes at the different 
offsets. In Figure 7 the top is the example embodiment, the middle is the NMO gather, 
and the bottom is a conventional migration technique. 

Figure 8 is a 2D model displaying the density. In this example embodiment, ray- 
tracing has been used to generate the synthetic data using velocities and densities. 
Figures 9 and 10 are the gather outputs of a conventional migration technique and an 
example embodiment of the present invention. The example embodiment of the present 
invention does a better job than the conventional technique in dip structure. 
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In a further embodiment of the present invention, a method is provided for 
migrating seismic data. The method comprises determining a travel time; determining a 
ray-bending correction; determining an amplitude correction; and outputting data based 
on the travel time, the ray-bending correction, and the amplitude correction. 

In a further embodiment, correcting amplitude further comprises restoring 
essentially relative true amplitude. In still a further embodiment, restoring relative true 
amplitude further comprises recovering angle dependent reflection coefficients. 

In a further embodiment, the method further comprises removal of a geometrical 
spreading factor. 

In an even further embodiment, determining an amplitude correction further 
comprises determining amplitude weights. In a further embodiment, the amplitude 
weights are dependent upon travel time. In a further embodiment, the amplitude weights 
further comprise emergence angles. In still a further embodiment the emergence angles 
further comprise: 



cos(o&) = 



1- 



f dT* 



\ 5X sJ 



V s 2 , and cos(tfg) = 



1- 



\ dx sj 



In an even further embodiment, the amplitude weights further comprise: 
Vcos^cos^g 



Vs 







T s ~ 


r 




Tl_ 



4L 2 H 2 sin* (a) 



+ 



2t/-4 



TT T'V 

s g r 



In an even further embodiment, the amplitude weights further comprise: 
Vcos^cos a s | det(rj Nsr + T^Nsr) \ 



v s V| det(JVs») |Vl det(iVc?*) | 
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In still a further embodiment, determining a travel time further comprises 
determining a velocity profile. 

In a further embodiment, determining a travel time further comprises determining 
an offset midpoint travel time. In a further embodiment, determining an offset midpoint 
travel time further comprises applying essentially the following formula: 



T.=k + X ** 



VL C,K 2 C+(1 + 2j7)j: ! ] 

wherein 

s-S 

where 8 and 5 are the anisotropic parameters. 

In a further embodiment, determining an offset midpoint travel time further 
comprises applying essentially the following formula: 



T = \ T « + yr~ c(VTi)xt + + -i - c(vti) x : 

V rms V ^ rms 

wherin 

CiVTI) = — 

In a further embodiment, determining a travel time further comprises determining 
a travel time for small offset to depth ratios. In a further embodiment, determining a 
travel time further comprises determining a travel time for large offset to depth ratios. 

In a further embodiment, determining a travel time for large offset to depth ratios 
further comprises applying essentially the following equation: 
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T=T +T +AT 

& s 



T 



2 T 2 



f i ^ 

\ + Lcc C4Xs 



2 T 



+ AT(VTI) 



In alternate embodiments, the method further comprises ray tracing and 
determining a travel time without ray tracing. 

In a further embodiment, the method further comprises correcting for VTI effects. 

In still a further embodiment, outputting the seismic data further comprises 
determining a time-table to calculate the output imaging time. In alternate embodiments, 
outputting the seismic data further comprises outputting the data in gather modes or in 
non-gather modes. 

In a further embodiment, determining a travel time further comprises determining 
a VTI travel time. In a further embodiment, determining a VTI travel time further 
comprises determining an offset-midpoint travel time. 

In a further embodiment, the determining and offset-midpoint travel time further 
comprises determining a migration velocity, and determining an anisotropy parameter. 

In a further embodiment, the method further comprises anti-alias filtering the data. 
In a further embodiment, the method further comprises determining a number of 
processors required to hold target bins. 

In a further embodiment, the correcting amplitude further comprises restoring 
relative true amplitude. 

In an even further embodiment, a system is provided for migrating seismic data. 
The system comprises means for determining a travel time. In alternate embodiments, 
means for determining a travel time comprise a computer, a workstation, any software 
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running on any computer, or any other means that will occur to those of ordinary skill in 
the art. The system further comprises means for determining a ray-bending correction. 
In alternate embodiments, means for determining a ray-bending correction comprise a 
computer, a workstation, any software running on any computer, or any other means that 
will occur to those of ordinary skill in the art. The system further comprises means for 
determining an amplitude correction. In alternate embodiments, means for determining 
an amplitude correction comprise a computer, a workstation, any software running on any 
computer, or any other means that will occur to those of ordinary skill in the art. The 
system further comprises means for outputting data based on the travel time, the ray- 
bending correction, and the amplitude correction. In alternate embodiments, means for 
outputting data comprise a computer, a workstation, any software running on any 
computer, or any other means that will occur to those of ordinary skill in the art. 
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In a further embodiment, the means for correcting amplitude further comprises 
means for restoring essentially relative true amplitude. In alternate embodiments, means 
for restoring essentially relative true amplitude comprise a computer, a workstation, any 
software running on any computer, or any other means that will occur to those of ordinary 
skill in the art. 

In a further embodiment, the means for restoring relative true amplitude further 
comprises means for recovering incidence angle dependent reflection coefficients. In 
alternate embodiments, means for recovering angle dependent reflection coefficients 
comprise a computer, a workstation, any software running on any computer, or any other 
means that will occur to those of ordinary skill in the art. 

In a further embodiment, the system further comprises means for removal of 
geometrical spreading factor. In alternate embodiments, means for removal of 
geometrical spreading factor comprise a computer, a workstation, any software running 
on any computer, or any other means that will occur to those of ordinary skill in the art. 

In a further embodiment, the means for determining an amplitude correction 
further comprises means for determining amplitude weights. In alternate embodiments, 
means for determining amplitude weights comprise a computer, a workstation, any 
software running on any computer, or any other means that will occur to those of ordinary 
skill in the art. 

In a further embodiment, the means for determining a travel time further 
comprises means for determining a velocity profile. In alternate embodiments, means for 
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determining a velocity profile comprise a computer, a workstation, any software running 
on any computer, or any other means that will occur to those of ordinary skill in the art. 

In a further embodiment, the means for determining a travel time further 
comprises means for determining an offset midpoint travel time. In alternate 
embodiments, means for determining an offset midpoint travel time comprise a computer, 
a workstation, any software running on any computer, or any other means that will occur 
to those of ordinary skill in the art. 

In a further embodiment, the system further comprises means for ray tracing. In 
alternate embodiments, means for ray tracing comprise a computer, a workstation, any 
software running on any computer, or any other means that will occur to those of ordinary 
skill in the art. 

In a further embodiment, the means for determining a travel time further 
comprises means for determining a travel time without ray tracing. In alternate 
embodiments, means for determining a travel time without ray tracing comprise a 
computer, a workstation, any software running on any computer, or any other means that 
will occur to those of ordinary skill in the art. 

In a further embodiment, the system further comprises means for correcting for 
VTI effects. In alternate embodiments, means for correcting VTI effects comprise a 
computer, a workstation, any software running on any computer, or any other means that 
will occur to those of ordinary skill in the art. 

In a further embodiment, the means for outputting the seismic data further 
comprises means for determining a time table to calculate the output imaging time. In 
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alternate embodiments, means for determining a time table comprise a computer, a 
workstation, any software running on any computer, or any other means that will occur to 
those of ordinary skill in the art. 

In a further embodiment, the means for outputting the seismic data further 
comprises means for outputting the data in gather modes. In alternate embodiments, 
means for outputting the data in gather modes comprise a computer, a workstation, any 
software running on any computer, or any other means that will occur to those of ordinary 
skill in the art. In a further embodiment, the means for outputting the seismic data further 
comprises means for outputting the data in non-gather modes. In alternate embodiments, 
means for outputting the data in non-gather modes comprise a computer, a workstation, 
any software running on any computer, or any other means that will occur to those of 
ordinary skill in the art. 

In a further embodiment, means for determining a travel time further comprises 
means for determining a VTI travel time. In alternate embodiments, means for 
determining a VTI travel time comprise a computer, a workstation, any software running 
on any computer, or any other means that will occur to those of ordinary skill in the art. 

In a further embodiment, the means for determining a VTI travel time further 
comprises means for determining an offset-midpoint travel time. In alternate 
embodiments, means for determining an offset-midpoint travel time comprise a 
computer, a workstation, any software running on any computer, or any other means that 
will occur to those of ordinary skill in the art. 
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In a further embodiment, the means for determining and offset-midpoint travel 
time further comprises means for determining a migration velocity; and means for 
determining an anisotropy parameter. In alternate embodiments, means for determining a 
migration velocity comprise a computer, a workstation, any software running on any 
computer, or any other means that will occur to those of ordinary skill in the art. In 
alternate embodiments, means for determining an anisotropy parameter comprise a 
computer, a workstation, any software running on any computer, or any other means that 
will occur to those of ordinary skill in the art. 

In a further embodiment, the system further comprises means for anti-alias 
filtering the data. In alternate embodiments, means for anti-alias filtering the data 
comprise a computer, a workstation, any software running on any computer, or any other 
means that will occur to those of ordinary skill in the art. 

In a further embodiment, the system further comprises means for determining a 
number of processors required to hold target bins. In alternate embodiments, means for 
determining the number of processor to hold target bins comprise a computer, a 
workstation, any software running on any computer, or any other means that will occur to 
those of ordinary skill in the art. 

In a further embodiment, the means for correcting amplitude further comprises 
means for restoring true amplitude. In alternate embodiments, means for restoring true 
amplitude comprise a computer, a workstation, any software running on any computer, or 
any other means that will occur to those of ordinary skill in the art. 



27 



ATTORNEY DOCKET NO. P30628US 

In an even further embodiment of the present invention, seismic data is provided. 
The seismic data is produced by the method comprising determining a travel time; 
determining a ray-bending correction; determining an amplitude correction; and 
outputting data based on the travel time, the ray-bending correction, and the amplitude 
correction. 

While many example embodiments of the present invention have been provided, the 
present description is presented to teach the broad aspect of the invention. This 
description is in no way exhaustive of the embodiments of the present invention. 
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